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We study extrema with velocity discontinuities of the variational two-body problem with bound- 
ary conditions in past and future. Along C 2 segments, these broken extrema satisfy Euler-Lagrange 
equations which are neutral differential delay equations with state-dependent deviating arguments. 
At points where velocities are discontinuous, broken extrema (spiky orbits) satisfy the Weierstrass- 
Erdmann corner conditions that energies and momenta are continuous. We construct periodic 
broken extrema near some C°° circular orbits, using solutions of the neutral differential delay equa- 
tions along regular segments and a variational approximation for the spiky boundary-layer segments. 
Finitely measured families of spiky periodic orbits result from the Weierstrass-Erdmann stabiliza- 
tion of the explosive tangent modes near some discrete circular orbits, with scales defined by the 
Weierstrass-Erdmann conditions. The Weierstrass-Erdmann conditions are further compatible with 
a condition to vanish the far-fields, thus creating a finitely measured family of periodic orbits with 
vanishing far-fields right at the atomic magnitude. We vary the mass ratio by three orders of 
magnitude, i.e., the hydrogen, muonium and positronium atoms. In each case the broken extrema 
have magnitudes in quantitative agreement with nature and quantum mechanics, i.e., the spacings 
between consecutive discrete orbital angular momenta agree with Planck's constant within thirty- 
percent, while orbital frequencies agree with a corresponding circular line of quantum mechanics 
within a few percent. The agreement is surprisingly better for the decay from the deepest quantum 
levels. Last, we discuss the extension of classical electrodynamics to variational electrodynamics. 

Keywords: calculus of variations, neutral differential delay equations, electrodynamics of point charges, 
state-dependent delay 



INTRODUCTION 

The variational formulation of classical electrodynam- 
ics [1] has continuous extrema with an undefined velocity 
at some points [2], henceforth broken extrema. Along 
C 2 segments, broken extrema satisfy neutral differen- 
tial delay equations with state- dependent deviating argu- 
ments [1, 2]. At corner points velocities/accelerations are 
not defined and broken extrema satisfy the Weierstrass- 
Erdmann corner conditions (WECC) that energies and 
momenta are continuous [3, 4]. We study the WECC for 
the variational electromagnetic two-body problem and its 
thresholds for corner creation. 

We construct periodic broken extrema of variational 
electrodynamics with a boundary-layer perturbation in- 
volving regular segments where the neutral differential 
delay equations and deviating arguments are linearized, 
while on the (thin) boundary-layer segments a variational 
approximation is used. Periodic broken extrema (spiky 
orbits), bifurcate from a discrete set of circular orbits 
with scales defined by the WECC. We vary the mass- 
ratio by three orders of magnitude to include the atoms 
of hydrogen, muonium and positronium. 

The Weierstrass-Erdmann mechanism of corner cre- 
ation stabilizes the tangent dynamics of spiky orbits, 
which are oscillations transversal to the particle sepa- 
ration. In the three cases of hydrogen, muonium and 
positronium, a discrete set of finitely measured orbits 



with vanishing far-fields have frequencies in agreement 
with quantum mechanics (QM) within a few percent. 
The qualitative agreement with QM is as follows: (i) 
the angular momenta of the broken minimizers are ap- 
proximately integer multiples of a basic angular momen- 
tum agreeing with Planck's constant within thirty per- 
cent, (ii) the emitted frequency is the difference of two 
eigenvalues of a suitable linear problem, just like in QM 
and (iii) the Weierstrass-Erdmann conditions involve pre- 
cisely the same relevant quantities of QM, i.e., continuity 
of momenta and energies at corner points. 

Maxwell's partial differential equations (PDE) repre- 
sent the experimentally verified integral laws of classical 
electrodynamics, e.g., Gauss's surface integral law for the 
electric field yields the divergence law when fields are C 1 . 
Differential laws were first written by Heaviside sometime 
after Maxwell completed electrodynamics in the 1860's, 
and good references exist [5]. Wheeler-Feynman elec- 
trodynamics of point charges (1945) [6] was constructed 
to avoid a shortcoming of Maxwell's PDE's, namely the 
divergence of the Lienard-Wiechert fields along charge's 
trajectories. Wheeler and Feynman [6] derived neutral 
differential delay equations of motion (NDDE), whose 
qualitative behavior was not understood at their time. As 
natural for the variational principle, solutions of NDDE 
must be defined piecewise [7, 8]. Continuation of solu- 
tions often leaves a set of zero measure where trajectories 
are not differentiable [7, 8]. 
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Quantum mechanics is an enormously successful prob- 
abilistic theory [9] , and since measure theory is the foun- 
dation of probability theory, the interface between vari- 
ational electrodynamics and QM is here consistently ex- 
plored with a measure-theoretic generalisation of classical 
electrodynamics . 

Electromagnetic fields of variational electrodynamics 
are not defined outside trajectories, and must be ex- 
tended with the Lienard-Wiechert formulas. Extrema 
with discontinuous velocities create discontinuous ex- 
tended fields, often undefined on sets of zero measure. 
Several results of Maxwell's electrodynamics do not carry 
over, e.g., Poynting's theorem is valid only in regions 
where the extended fields are C 2 . Gauss's law can be 
shown by density arguments to hold for extended fields 
of trajectories having a finite number of corners [2]. Ex- 
tended fields of trajectories with an infinite number of 
corners per average time unit have not yet been studied. 
Generalizations of the integral laws of classical electrody- 
namics using Sobolev's trace theorems[10], and measure- 
theoretic generalizations to Lebesgue-integrable action 
functionals, are also open problems. 

The electromagnetic two-body problem has a known 
one-parameter family of C°° solutions, i.e., the circular 
orbits [11], which create non- vanishing far-fields and are 
minimizers of variational electrodynamics as well. A sur- 
prising result for two-body orbits with repulsive inter- 
action and motion restricted to a line [12] is that glob- 
ally defined C 2 orbits are uniquely determined by the 
distance of closest approximation, i.e., a one-parameter 
invariant manifold. The analogous problem within vari- 
ational electrodynamics has not been studied yet. The 
stability of circular orbits was first studied in [13, 14] and 
[15]. Action-at-a-distance is discussed in [1, 16-18], some 
electromagnetic-like delay equations are discussed in [19] 
and the WKB quantization of circular orbits is discussed 
in [20, 21]. 

This paper is divided as follows: In Section II we 
explain variational electrodynamics as an extension of 
Wheeler-Feynman electrodynamics. We introduce the 
variational structure, the WECC, and explain the gen- 
eralized absorber condition needed in variational electro- 
dynamics. In Section III we introduce the Schild circular 
orbits and magnitudes for later use, which is our starting 
point to construct a fast-scale (stiff) tangent dynamics. 
In Section IV we linearize the Wheeler-Feynman NDDE 
about circular orbits. We explain in detail the infinite- 
dimensional tangent dynamics of circular orbits and the 
infinite number of linearly unstable transversal modes. 
Section V contains the boundary layer theory and the 
Weierstrass-Erdmann threshold for corner creation. In 
Section VI we study the magnitudes prescribed by the 
WECC for a finitely measured family of periodic broken 
extrema. In Section VI we compare the magnitudes of 
the broken extrema with nature and QM for hydrogen, 
muonium and positronium. In Section VII we discuss 



the GAC mechanism to vanish the far-fields and show 
how it is surprisingly compatible with the WECC. Last, 
Section VIII contains discussions, speculations and the 
conclusion. 



VARIATIONAL ELECTRODYNAMICS 

Here we discuss an extension of Wheeler-Feynman's 
classical electrodynamics of point charges [6]. Unlike 
Maxwell's electrodynamics, variational electrodynamics 
does not involve the solution of a PDE: Instead, finite 
segments of charge's trajectories are critical points of a 
variational principle with Einstein-local boundary condi- 
tions in past and future [1, 2}. AC 2 extremum of varia- 
tional electrodynamics satisfies Euler-Lagrange equations 
where the Lienard-Wiechert fields of other charges en- 
ter as coupling terms. These Euler-Lagrange equations 
hold only along the finite segment of trajectories of the 
variational problem. Elsewhere in space electromagnetic 
fields are obtained by extension, often modulo a set of 
zero measure in light-cone with breaking points. Exten- 
sion of fields to almost everywhere in R x M 3 requires 
globally defined extrema, otherwise extension is limited 
inside a ball of finite radius. 

We henceforth use a unit system where the speed of 
light is c = 1 and the electronic and protonic masses and 
charges are mi = 1 and qi = —1 and mi — 1836.1526 and 
q2 = 1, respectively. Electrodynamics until the 1900's 
had in mind the Lorentz force equation of motion for 
charge qi , i.e., 



dt X^/T^ 2 



ft(E(x f ,t)+v i xB(x jl t)) 1 (1) 



where particle i has mass and x^(i) e C 2 (M.). The 
idea of solving Maxwell's PDE's for the fields E(xj, t) and 
B(x,j,i) entered into a crisis when the Lienard-Wiechert 
general solution for the field of a system of point charges 
came out [5], i.e., 



E(xi,t) = ^(aE i+ + 6E J -_), 



j 



B(xi,i) = > (aB j+ + 6B,_), 
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(2) 
(3) 



where summation is over all charges, a and b are any pair 
of real constants satisfying a + b = 1 [5] and 

F - n < Uj± I nj± X (Uj± X &J±) l (A\ 



j± ~ (l±n j± .v i± )3' 
B J± = Tn J± x E j± , 

-1/2 



(5) 
(6) 



with7j± = (l-v 2 ± ) 1/2 , Vj± = dx.j/dt\ t =t j± and a 3± = 
d 2 Xj /dt 2 \t=t j± the velocity and acceleration of charge j 
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evaluated at the retarded/advanced time tj^ denned by 

*iT=*T||x i (t J - T )-x)|| J (7) 

an implicit state-dependency on x.j(tj^) henceforth 
called the light-cone condition. Last, in Eqs. (4), 
rij± = (x — ~x-j(tj^))/ri± is a unit vector from position 
Xj(tj^) to point x and rj± = ||x 



JT) 



denotes dis- 
tance in light-cone. Maxwell's PDE's fields (4) diverge 
along each charge's trajectory: Wheeler and Feynman's 
recipe for a finite Lorentz force [6] was to restrict the 
summation in Eqs. (2) and (3) to j =/= i, excluding the 
divergent self-field. 

The choice i 7^ j combined with a = 6=l/2is the 
solution of the variational boundary value problem ex- 
plained bellow. Even restricted to almost everywhere 
subluminal, non-collisional trajectories [1], the electro- 
magnetic action principle is more complex than the two- 
point boundary value principle of Newtonian mechanics 
since boundary segments in past and future must be given 
to provide data for the advanced and retarded arguments 
[1], as illustrated in FIG.(l). The needed Einstein-local 
boundary conditions consist of (a) the initial point Ox for 
trajectory 1 and respective boundary-segment of trajec- 
tory 2 inside the light-cone of 0\ , and (b) the final point 
L 2 for trajectory 2 and respective boundary-segment of 
trajectory 1 inside the light-cone of L 2 . The variational 
problem is to find the segments of trajectories between 
the endpoints Ox and L 2 , and Of and L 2 of FIG. (1). As 
shown in [18] and illustrated in FIG. (1), along almost 
everywhere subluminal continuous orbits satisfying these 
boundaries, the light-cone conditions (7) define unique 
solutions for the deviating arguments tj±(t, Xj(t)). 

The action functional is a sum of four integrals: 
two local integrals each involving one trajectory, i.e., 
J Ti(xj,Xj)df,*, and two interaction integrals of func- 
tions of both positions/velocities, where one posi- 
tion/velocity is evaluated at a deviating time argument, 
/Vy (xj,Xi,Xj±,Xj±)d£j. The action can be expressed 
in two equivalent forms, i.e., 



S = [ L2 T 2 dt 2 + f ^Txdh 
Jt n + Jt 0l 



2 , -r — 



V-^dtx+J V^dh, 
to, Jto 1 



(8) 



Txdtx 



to 1 



T 2 dt 2 



V+dt 2 



V 21 dt 2 , 



(9) 



where vertical braces indicate equivalence by change of 
the integration variable on the interaction integrals. The 
Jacobian for the change of variables is the derivative of 
the delayed times tj±(ti,Xj(ti)) evaluated along the or- 
bit, as obtained taking a derivative of the implicit condi- 




FIG. 1. Schematic showing the Einstein-local boundaries in 
R 4 , i.e., (a) initial point Oi = (to 1 ,'Xi(to 1 )) of trajectory 1 
and boundary segment X2^) for £2 G [t -,t Q +] inside the 
light-cone of Ox, and (b) endpoint L2 = (ti 2 , X2(£l 2 )) of tra- 
jectory 2 and boundary segment Xi(ti) for tx £ \t L - , t L +] in- 
side the light-cone of 1/2. Trajectories xi(ti) for tx G [<Oi >"t L -\ 
(green line) and 'x.%{t%) for ti 6 [t + , tz, 2 ] (blue line) are deter- 
mined by the extremum condition. Black dotted lines indicate 
the light-cone relation and brown dotted lines indicate the 
and interactions. Arbitrary units. 



tion (7) with x = Xj(tj), 



dtj± _ (1 ± n J± ■ Vj) 



dt; 



(1 ± n J± ■ vj±) 



(10) 



and explained in Refs. [1, 18]. One can thus express the 
interaction terms by either integrals over tx (as in (8)) or 
integrals over t 2 (as in (9)). A generic variational struc- 
ture can be defined using arbitrary V's on line (8), while 
the V's on line (9) are determined by changing variables 
with (10) or the equivalent Jacobian if the constraint is 
other than (7). Henceforth we concentrate on variational 
electrodynamics, i.e., 



1 



V?. 



(l-Vi-V,-±) 



2r i± (l±n i± ■ v 3 -±)' 



(11) 
(12) 



where i = 1,2 and j = 3 — i. 

For the linear variation, trajectory 2 is to be kept fixed 
while trajectory 1 is varied, and line (8) with the first 
term constant defines partial Lagrangian 1. Vice- versa, 
trajectory 1 is to be kept fixed while trajectory 2 is var- 
ied, and line (9) with the first term frozen defines par- 
tial Lagrangian 2. AC 2 extremum satisfies the Eulcr- 
Lagrange equation of each partial Lagrangian [1], i.e., 
Eq.(l), the NDDE with state-dependent deviating argu- 
ments of Wheeler and Feynman[6]. 

Next in complexity are broken extrema with a finite 
number of corner points [4] . For each partial Lagrangian, 
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broken extrema satisfy the WECC that momenta and 
energies are continuous at corner points. The first 
Weierstrass-Erdmann corner condition, (WECC1), is the 
continuity of the momentum of partial Lagrangian i, i.e., 



<9v,; 



v/T^? 2r i _(l-D i _.v J -_) 
2r J -+(l + n i+ - V j+)' 

2 J' 



(13) 
(14) 



where Jz?i is the sum of the last three integrands of line 
(8), Jt?2 the sum of the last three integrands of line (9) 
and Aj the vector potential of particle j. Equation (13) 
holds for i = 1,2 despite any velocity discontinuities of 
charge j = 3— i. If the extremum has a corner for particle 
i at time £j, the trajectory of particle j compensates with 
a corner in light-cone at either tj- or tj+. 

For the second Weierstrass-Erdmann corner condition, 
(WECC2), we express the partial Lagrangian i as 



(15) 



where 
A., E 



*3 + 



2r i _(l - nj_ • Vj-_) 2r J+ (l + n j+ ■ v j+ ) ' 
1 1 



3 2rj-(l - nj_ • Vj -_) ' 2r J+ (l + n j+ ■ v j+ ) ' 
The partial energy of partial Lagrangian i evaluates to 



E, 



<9v, 



(16) 



Therefore, the WECC2 is simply that |vj| must be con- 
tinuous at each breaking point of each trajectory, in order 
to make the partial energies Ei continuous. 

Next we express the WECC1 using an upper index I 
or r to indicate respectively left- velocity or right-velocity 
at the breaking point. Using (13) and the WECC2, Eq. 
(16), the WECC1 becomes 



Avj_ - rr, x (v£_ x v^-_) 
;l<\, ;1 n, -v- iil n, -v; 
Av j+ + n i+ x ( V J + x vj + ) 



^2r J+ (l + n J+ -v' + )(l + n J+ -v J r 



= 0, (17) 



where Avi = v[ - v|, Av^± eee vj ± - vj ± . 

Last, we already mentioned that fields of trajectories 
defined in t G [—00,00], henceforth globally bounded, 
can be extended to almost everywhere in R x M 3 . By 
extension we mean fields defined by formulas (4) and 



(6) whenever all the respective retarded/advanced veloc- 
ities/accelerations are well defined. At least for trajecto- 
ries with a finite number of corners per unit time, formu- 
las (4) and (6) define extended fields almost everywhere 
but on a set of zero measure. Extension is used when an- 
other charge is added to the system, at which stage min- 
imization respect to its trajectory yields an extra Eulcr- 
Lagrange equation involving the extended fields. Let 
us first assume globally bounded trajectories £ C 2 (R), 
which implies that all advanced and retarded normals 
coincide at a far distance, i.e., n J+ = rij- = n, V j [2], 
yielding the extended electric field 



E(x,t) 



(x,t) + E,_(x,t)), (18) 



while Eq. (6) with rij + = rij = n yields 

B(x,t) = \n x ^(E 3 -_(x,t) - E j+ (x,i)). (19) 

j 

According to Eq. (1), for an extra distant charge to travel 
undisturbed at a far distance with an arbitrary velocity 
Vi, the semi-sum (18) and semi- difference (19) of univer- 
sal electric fields must vanish asymptotically. Wheeler 
and Feynman's absorber condition was obtained using 
Maxwell's PDE's [6, 22], while the above pre-minimizer 
condition applies also to piecewise C 2 globally bounded 
extrema, yielding the separate vanishing almost every- 
where of the retarded and advanced universal far fields. 
We henceforth call this the generalized absorber condition 
(GAC) of variational electrodynamics. 

Notice that functional (8) is semi-bounded (S > 0), 
and under some compacity properties one extremum 
should be a minimum [23]. It is an open question if 
(8) could also have saddle-point type extrema [18]. 

In physics, the name classical trajectory is synonymous 
of global solution of an ordinary differential equation 
(ODE), i.e., a set of points invariant by time evolution 
under the ODE. In variational electrodynamics, the anal- 
ogous concept is a global C° curve out of which one can 
arbitrarily cut finite boundary segments, as defined above 
FIG. 1. The invariant manifold can have corners and it 
is only necessary that one extremum of each finite vari- 
ational problem be a segment of the manifold. Unlike 
the case of ODE's, one can not continue an invariant 
manifold forwards (nor backwards) with a time integra- 
tion. Instead, the infinite dimensional problem at hand 
is finding a whole function, just like solving a PDE with 
boundary conditions. 

Last and again, trajectories with derivative disconti- 
nuities obviously do not have a Taylor series approxima- 
tion. Together with the failure of Poynting's theorem 
for extended fields of spiky orbits, these are two ma- 
jor generalisations of Dirac's classical electrodynamics of 
point charges [24]. The generalization is welcome after 
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the proof of Ref. [2] that globally bounded two-body 
orbits with vanishing far-fields must have corners. 



CIRCULAR ORBITS 

The details of the C°° Schild circular orbits [11] are 
given below using the notation of Ref. [15]: The constant 
angular velocity and distance in light-cone are denoted 
by Q and r&, respectively, and the angle 9 = firt, that 
particles turn in the light-cone time is henceforth called 
the delay angle. The family of subluminal C°° circular 
orbits [11] is parametrized by 9 G (0, 1), and for atomic 
orbits 9 < 10~ 2 . Along orbits with a small delay angle 
9, the Kepler formulas yield the leading order angular 
velocity and distance in light-cone, respectively, i.e., 



Q = fi9 3 + 0(9 5 ) 
9 1 



n 



n fie 2 ' 



(20) 
(21) 



where reduced mass and total mass are defined by fi = 
m\mijM and M = mi + m^. It is important to keep 
these dependencies in mind, and that for hydrogen fi ~ 
(1836/1837) ~ 1. Adopting the same notation of Ref. 
[15], we express each particle's circular orbit's radius by 
numbers < 6, < 1 as 




FIG. 2. Schild circular orbit with particles in diametral oppo- 
sition at the same time of the inertial frame (dotted straight 
line is t = 0). Indicated is the retarded position of particle 
1 along the outer circle, and the angle 8 travelled during the 
light-cone time. Particle separation (rj) in the past light-cone 
of particle 2 is the red line making an angle ip2 with the dotted 
line. The delay angle and the protonic circular radius have 
an exaggerated magnitude for illustrative purposes. Arbitrary 
units. 



LINEAR MODES 



which define scalar velocities 

v t = ftr l = 9bi 



(22) 



(23) 



for % = 1,2. 

As illustrated in FIG. (2) and in Ref. [11], the circular 
radii and the distance in light-cone form a triangle of 
largest side r^, yielding a trigonometric constraint 



26i6 2 cos(6») = 1, 



(24) 



equivalent to Eq. (3.1) of Ref. [11]. 
of equations (3.2) and (3.3) of Ref. 



Using the ratio 
[11] to eliminate 



O, together with constraint (24), we find the solution 



' 2M l 



= (1 + ^)^+0(0 and 6 2 = (1 
Last, Noether's angular momentum 



1 + b x b 2 9 2 cos(0) _ 1 
9 + b 1 b 2 9 2 sin(0) ~ 9' 



0(9 4 



(25) 



is an important quantity of the circular orbit to keep in 
mind [11, 14, 20]. Atomic orbits have l z of the order of 
one over the fine-structure constant, hc/e 2 = 137.036, 
a fundamental magnitude of nature. Delay angles and 
angular momenta of Bohr orbits in our unit system are 
respectively 9 ~ (e 2 /hc)/q and l z ~ (hc/e 2 )q, for each 
nonzero integer q. 



In order to linearize the equations of motion about cir- 
cular orbits of small delay angles, we re- write the Lorcntz- 
force Eqs. of motion (1) as 



qi[F,j - (Vj • Ej)v 4 + v f x Bj], 



(26) 



by evaluating the derivative on the left-hand side of Eq. 
(1) and subtracting the scalar product with Vj, where i = 
1, 2 and j = 3 — i. The magnetic term (last term on the 
right-hand side of (26)) is a transversal force proportional 
to the electric field by (6) , and further proportional to the 
velocity, which is small for small 9 (see Eq. (23)). 

The electric force, first term on the right-hand side of 
Eq. (26), is responsible for the significant contributions 
to the linearized equations along a circular orbit: The 
electric field (4) decomposes in two terms: (i) the near- 
electric field proportional to V r f±) a magnitude of /i 2 9 4 
by use of (21), times an (almost) radial vector, and (ii) 
the far-electric field proportional to l/rj±, a magnitude 
of [i9 2 by use of (21), times a vector (almost) along the 
circular orbit. It is important to ponder upon the magni- 
tude of each electric contribution: the far-electric field is 
(almost) transversal, and projection to the normal direc- 
tion brings another factor of 9. The force along circular 
orbits is approximately described by the near field only 
because the far-field (4) is further proportional to the ac- 
celeration (which is 0(l/r 2 ) along circular orbits). The 
contributions to the right-hand-side of (26) have magni- 



G 



tudes 



Fi oc 0(6 4 ) + |a J -|0(6» 3 ). 



(27) 



In Eq. (27), the near-field contribution is 0(9 4 ) while the 
far-held contribution is \&j\0(9 3 ), which is smaller along 
circular orbits since |aj| oc 1/r 2 = 0(9 4 ). Upon lineariza- 
tion about the circular orbit, this dominance changes: the 
linearized equations accept solutions of arbitrarily large 
frequencies, and the most important linearized term is 
precisely the far-field term. 

Next we derive the linearized equations along the or- 
bital plane using the notation of Ref. [15]. We introduce 
complex gyroscopic coordinates where the circular orbit 
is a fixed point of the equations of motion, i.e., 



+ Wj = n exp(- 



-iClt)[bj + lj 



(28) 



where (lj,Uj) are respectively the longitudinal and 
transversal gyroscopic coordinates. The circular orbit 
[11] is the fixed point (lj,Uj) = (0,0) for j = 1,2. The 
dominant linear correction for large accelerations (the 
stiff limit [15]) is obtained from (26) with only the first 
term on the right-hand side and far-field Ej± approxi- 
mated by 



E 



3± 



— n i± x(n,-±x<yaj±) = —{n d ±-6a J ±)n j ±-—5aj±, 

' b 1 b 1 b 

. (29) 

where ± indicate evaluation at the unperturbed deviating 
times t ± (O/Cl) (because we want the linear term only). 
The gyroscopic representation of the rotating normals 
in light-cone are complex numbers of unit modulus, i.e., 
exp(— iilt) exp(iipj), where 



exp(i^i) = 61+62 exp i6, 
exp(itp2) = 62 + ^1 exjpiO. 



(30) 



The unitarity of (30) is seen using Eq. (24). Angles ip\ 
and ip2 further satisfy ipi + ip2 — $7 and for small 9 the 6, 



given below Eq. (24) yield ip% 



M 



9 and i/j 2 



M 



In 



FIG. (2), ip2 is the angle between the (rotating) red line 
and the a;- axis (dashed line). 

The linearized planar equations of motion are obtained 
substituting (28) into (26) with only its first right-hand- 
side term given by the semi-sum of retarded and ad- 
vanced fields (29). The real and complex parts of the 
linearized equations of motion keeping only the largest 
derivatives of the gyroscopic coordinates yields 



m 2 r b i 2 
m 2 r b u 2 



(I2+ + I2-) G (&2+-U2-) 
--Cu g 611 2 ' 

(h+ + 'h-) (tii+-iii_) 



-C 



12" 



S12 ' 



(U2+ + U2-) „ {h+-h-) 

-C 13 - 6 13 - , 

(ui+ + ui_) 



-Ci 



Si 



(31) 



where C±j and S%j are given below Eq. (33). Equation 
(31) is a linear NDDE with exponential solutions 

l = Ljexp(\tit/9), 
uj = UjeyLp(\nt/9), (32) 

1,2 and L\, L2, U%, U2 are non-trivial solu- 



where i, j 
tions of 



mxrt, Dn En 

D12 m 2 r b E12 

E13 m 1 r b D 13 
En D14 m 2 r b 




where 



(cos Q — cos ip2 cos £ 
cos ^+sin ipi sin 6 



and 



(sin 9— sin tp2 cos 9 s 
sin 6— sin thi cos 9 
(l-cosV-2)sine 
(1— cos tpi) sin 9 t 



(33) 



(34) 



(35) 



A nontrivial solution of (33) requires the vanishing of 
its 4 x 4 determinant, 



F xv = l- ^cosh 2 (A)=0, 



(36) 



where powers of 9 2 cosh A with coefficients of 0(9 4 ) were 
discarded, henceforth called the stiff-limit [15]. The 
roots of Eq. (36) exist in symplectic sets of four, i.e., 
(A, -A, A*, —A*). For atomic hydrogen, - 137.036 
and (jji/M) ~ (1/1836), so that Eq. (36) requires A to 

have a real part |5R(A)| = a = \n{J^) ~ 14.29, while 

the imaginary part of A is an integer multiple of iri. The 
general solution of Eq. (36) modulo the symplectic sym- 
metry is 



A = a 



nqi, 



(37) 



where i = y/—l and q E Z. Solution (37) was called a 
ping-pong mode in Ref. [15] because its phase advances 
by irq in one light-cone time r b , i.e., a phase speed of 
nqn/9 = n9 2 irq = nq/r b = (f7/6»)3(A). The only 0(1) 
off-diagonal terms of matrix (33) are D13 and D14, others 
being 0(9) or higher order. The nontrivial eigenvector 
solution is approximately 



(L 1 ,L 2 ,C7 1 ,f/ 2 )cx( 



fi6 n6 



(38) 



so that for m 2 3> m\ solution (32) is along the circu- 
lar orbit, since Ui 3> Lj, i.e., a quasi- transversal mode. 
The largest longitudinal component is for positronium 
at moderate 9, but still small. Solution (37) defines a 
nonzero real part for A, so that amplitudes blow up at ei- 
ther t — > ±00, implying that besides circular orbits there 
are no other C°° almost-circular orbits defined for all 
times (i.e., globally bounded). 
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The inclusion of 0(1/ X) and 0(1/ X 2 ) linear terms to 
the planar motion is outlined in Ref. [15]. The lin- 
earized motion perpendicular to the orbital plane is stud- 
ied analogously. As explained in [15], the z-oscillations 
are transversal modes that decouple from the planar 
transversal oscillations (31) at linear order. The de- 
terminant of the linearized 2x2 system is again (36), 
an asymptotic degeneracy. The degeneracy is raised by 
0(1/X 2 ) corrections to (33) and (36) introduced by the 
linear terms with lower derivatives [15]. The determinant 
for planar modes is calculated in Ref. [15] up to O(i) 
terms (see Eq. (41) of Ref. [15] with V = -1/2 ), i.e., 



F. 



1 



f^ 4 



(1 



7 



A2 +^) cosh 2 (A) 



iif) 4 1 5 
+ ^ ( A + A^ )sinh(2A) = ' 



(39) 



while the determinant for perpendicular oscillations and 
up to O(ji) terms is Eq. (B17) in appendix B of [15] 
with T = -1/2, i.e., 



1 



(1 



M { X 



1 

~ A2 



1 

+ F 
sinh(2A) 



) cosh 2 (A) 

0. (40) 




FIG. 3. Unperturbed circular orbits (dashed lines), boundary 
layers along orbit of (heavier) particle 1 (solid red) diving to a 
corner at a smaller inner radius (dotted line). Boundary layers 
along orbit 2 (solid blue) spring up to corners at a larger outer 
radius (dotted line). At points (h), (a) and (c), perturbed 
orbit 1 (dark solid line) crosses the unperturbed orbit and 
makes a regular segment with an exaggerated outer radius. At 
points (d) and (f), perturbed orbit 2 (inner dark circular line) 
crosses the unperturbed orbit and makes a regular segment at 
a smaller radius. Boundary layer magnitudes are exaggerated 
for illustrative purposes. Arbitrary units. 



Notice the symplectic symmetry that roots of Eqs. (39) 
and (40) are still in sets of four, (A, —A, A*, —A*), and 
even the first correction at O(j) is the same. The cor- 
rections at O(jz) separate Eq. (39) from (40), unfolding 
the asymptotic degeneracy of planar and perpendicular 
transversal modes. 



BOUNDARY LAYER 

The non-zero real part in Eq. (37) suggests that there 
are no other near-circular C°° solutions besides circu- 
lar orbits, since all linearized modes blow up at either 
t — > ±oo. This unless a corner point invalidates the lin- 
earization. Moreover, Ref. [2] already found that globally 
bounded two-body orbits with vanishing far-fields (the 
GAC) must have velocity discontinuities. Therefore, we 
are led to seek globally bounded extrema with corners. 

We set about the task to construct a periodic bro- 
ken extremum using a boundary-layer perturbation that 
assumes a number of piecewise C°° layers separated 
by boundary layers each containing one or more corner 
points, as illustrated in FIG. (3). Along both trajectories, 
boundary layers have an angular width aO <§C 9. Outside 
boundary layers we can linearize because trajectories are 
C°° and deviating arguments fall in C°° segments. In- 
side boundary layers we do not linearize but rather use a 
variational approximation for the minimizer with corners 
and apply the WECC. 



The exact roots of Eqs. (39) and (40) near the limiting 
root (37) are, respectively, 

X xy (0, q) = o~ X y + nqi + ie xy , (41) 
X z (0,q) = a z + irqi + ie z , (42) 

where q is an arbitrary integer, t xy (9,q) = —e xy (9,—q) 
and e z (9,q) = —c z (9 1 —q) (by the symplectic symmetry). 
An 0(/j,9 4 /M) approximation for positive a is obtained 
by setting 4 cosh 2 (A) ~ 2sinh(2A) ps exp(2A) into Eqs. 
(39) and (40), thus defining polynomials 

7 10 5 

A2 ~ A3 + W 



F X y(X) 



ri9 4 , , w 2 



and 



F Z (X) = 1 



= 1 



^ex P (2A)(l---- 
^ex P (2A) Pz (-) = 0. 



2 

A^ 



It follows from Eqs. (39) and (41) that 

Pxy ( Jir- ) 



exp(4ie xy ) 
while Eqs. (40) and (42) yield 



exp(4ie z ) 



*(£) 



(43) 



+ A^ 



(44) 



(45) 



(46) 
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Disregarding the contribution of e xy and e z and replac- 
ing X xy = X z = a + nqi in the right-hand-sides of (45) 
and (46) yields, to the first order in (1/A) 



-irq 



(a 2 + w 2 q 2 ) 



(47) 



We henceforth use a positive integer q, so that the ener- 
getic mismatches e xy and e 2 are negative O(j) numbers. 
Another order can be gained expanding the right-hand- 
side of (45) in a Taylor series on the deviation e xy about 
\ q = a + nqi, which generates only terms at 0(1/A 3 ), so 
that up to 0(1/A 2 ) Eq. (45) yields 



(?) 



-nq 



bnqa 



{a 2 +ir 2 q 2 ) (a 2 + ir 2 q 2 ) 2 ' 



(48) 



e xy negative and monotonically increasing for q > 6. 
Analogously, the right-hand-side of (46) evaluated at 
X q = a + nqi yields e z up to 0(1/A 2 ) 



M = 



-nq 



Sirqa 



(a 2 



n 2 q 2 ) 



{a 2 + 7T 2 q 2 ) 2 ' 



(49) 



again negative and monotonically increasing for q > 6. 

Let us assume a corner at t = along orbit 1, as il- 
lustrated in FIG. (3). We define the first regular layer 
art, < t < r& — art, as the first and last zeros of the (ex- 
ponentially increasing) perturbation inside the first light- 
cone time zone. This is illustrated in FIG. (3) as points 
(a) and (h) where the perturbed trajectories (solid lines) 
cross the unperturbed orbits (dotted lines). From (a) to 
(h) the velocity of the perturbed orbit changes almost by 
irq. Using a linear combination of the symplectic quartet 
of linearized modes, and (28) and (38) with li cx L\ = 
and u\ cx U\ = 1, the orbital perturbation constructed 
to vanish at layer edges is 



(f) 



cosh(^) 
cosh(^p) 



sm 



[irq + e xy (q) - 9](t - ar b ) 



n 



(50) 



In Eq. (50), a xy and e xy are given by (41) and we must 
have the amplitude A < to avoid a superluminal 
velocity at layer edges. Along the regular region the 
phase of the sine function in (50) advances by almost 
nq, and the condition of vanishing at t = r& — ari, yields 
a - 

An instructive variational approximation for the 
boundary layer motion is illustrated in FIG. (3), with 
corners along a given orbit all equivalent by a ^-rotation. 
The central spike starts along trajectory 1 with a ninety- 
degree corner, point (a), and respective corner in light- 
cone along trajectory 2, point (d). Next is a straight-line 
dive (ab) to an almost 180-degree corner ((b) and (e)), 
then a straight-line climb (be) to the last ninety-degree 



corner ((c) and (f)), resuming motion along the next reg- 
ular layer. 

Corners satisfy the WECC, and 180-degree corners, 
(points (a) and (d) in FIG. (3)), are simpler to analyze: 
For planar motion, condition (17) yields four equations, 
one along each Cartesian direction and for i — 1,2. As- 
suming 180-degree corners are reflections about the light- 
cone direction n, the vectorial components of (17) per- 
pendicular to each n vanish trivially by symmetry. The 
two nontrivial conditions yield a 2 x 2 linear homogeneous 
system for the n • Avj , with determinant 



D = m\mi\ 1 



1 



(51) 



where f is distance in light-cone from (a) to (d). In de- 
riving (51) from (17) we assumed 180-degrees reflections, 
otherwise there is a contribution from the (vj ± x vj ± ) 
terms. Ninety degree corners are analyzed substituting 
v l = v \ + m to (17), yielding a 4 x 4 linear system 
for the Avj, of which (51) is a special case. 

Equation (51) poses an alternative: if distance r is to 
be near the circular light-cone distance = -p-, then 
the last term on the right-hand-side of (51) is negligi- 
ble and the vanishing of the determinant requires that 
at least one particle has a quasi-luminal velocity. Oth- 
erwise collisions can happen only at distances near the 
very small protonic radius. In variational electrodynam- 
ics, 180-degree corners function as stepping-stones when 
linearizations reach the critical quasi-luminal velocities 
necessary for the formation of corners at-a-far-distance. 
Above threshold, corner creation is a rebouncing mech- 
anism that offers an alternative to charges falling into 
each other. The stepping-stone condition further implies 
that the amplitude of transversal mode (50) at layer edge 
must be A ~ nL . 

irq 

Recalling that motion in the regular layer is almost- 
transversal, the WECC2 (16) poses another threshold 
for the 90-degree corner followed by a straight-line quasi- 
luminal dive (point (a) in FIG. (3)): since | Vx | must be 
continuous, the dive requires a quasi-luminal transversal 
velocity at layer edge by continuity. Again, the exponen- 
tial blow-up of Eq. (50) reaches a quasi-luminal velocity 
at layer edge if A ~ ^ . 

FIG. (4) illustrates the boundary layer motion in the 
presence of planar and perpendicular perturbations. The 
green lines of FIG. (4) indicate corners connected by the 
light-cone condition, again with exaggerated magnitudes 
to illustrate the gapped-solenoidal motion. The layering 
must fit an integer number of regular regions and bound- 
ary layers in the circle. Assuming the time of flight for 
each abc segment of FIG. (3) to be At ~ 2rf,, the total 
time for N c corners is 2a — 2iV c r;,, yielding the rough 

estimates Q ~ = {8 ~^ y) . 
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FIG. 4. Boundary layer motion along a 3-dimensional orbit 
with perpendicular and transversal perturbations. Perturbed 
trajectory of (heavier) particle 1 (solid red) and perturbed 
trajectory 2 (solid blue). Solid green lines indicate corners in 
light-cone and dashed dark lines are the unperturbed orbits. 
Boundary layer and delay angle have exaggerated magnitudes 
for illustrative purposes. Arbitrary units. 



WEIERSTRASS-ERDMANN AND 
PROBABILITY 

The tangent dynamics of circular orbits has an infi- 
nite number of unstable transversal modes of arbitrarily 
large frequencies, as seen by the non-zero real parts of 
linear growths (41) and (42). This is unlike the classical 
Coulombian ODE approximation for the two-body prob- 
lem, whose tangent dynamics has a finite number of fre- 
quencies of the order of the orbital frequency [25]. Since 
linearized modes (41) and (42) are unstable, continuation 
along the C°° segment would blow up, and a velocity dis- 
continuity is needed to break from the C°° segment. A 
corner requires stepping-stone condition (51), and for a 
finitely measured family of orbits with corners to exist, 
the WECC2 along orbits with planar and perpendicular 
perturbations requires the resonances studied below. 

The perpendicular perturbation in the regular region is 
constructed analogously to the planar perturbation (50), 
using the linearized modes explained in appendix B of 
Ref. [15], and yielding 



cosh(— t) 
cosh(^) 



[its + e z (s)](t - ar b ) 



, (52) 



where s is any integer (possibly different from q) and B is 
an arbitrary amplitude smaller than ^ to avoid a super- 
luminal velocity at layer edge. Again, the amplitude of 
the transversal perturbation (52) is constructed to van- 
ish along the circular orbit at t — art and t = rt — art, 
to keep the property that corners see other corners in 
light-cone (a boundary-layer-adjusted resonance). 

Since the phase of oscillation (~ ns/r b ) is fast, a large 
orbital velocity at layer edge results from a small ampli- 
tude B ~ which is illustrated in FIG. (3). Layer edge 



amplitudes of both types of transversal modes should 
vanish while their derivatives reach a quasi-luminal ve- 
locity. Amplitude (50) vanishes at t = r& — art, when 



[nq + txy(q) -9](1- 2a) = irq, 



(53) 



other multiples of ir being excluded because e xy (q), 9 
and a are small. Analogously, amplitude (52) vanishes 
at t = r\, — art when 



[ivs + e z (s)](l - 2a) = ns, 



(54) 



where again other integer multiples of 7r arc impossible 
because e z (s) and a are small. 

The climbing/diving edges at t = art, and t = Vb — art 
of the 'periodic broken minimizer are equivalent, imply- 
ing that (50) and (52) with arbitrary amplitudes near 
(A,B) — (-^S J^) should yield the same quasi-luminal 
|v., |. For this it is necessary that 



qexy(q) - se z (s) 



(55) 



Condition (55) ensures the WECC2 that \vj\ is contin- 
uous at corners of the periodic broken minimizer and is 
a probabilistic argument. Given (55), at points (a), (c) 
and (h) of FIG. (3) a finitely measured family of orbits 
include a corner turn, i.e., a whole neighbourhood of per- 
turbed orbits near (A.B) = ( — , —). Next we calculate 
the magnitudes of these resonant orbits with q = s. 

For each q = s, condition (55) determines a unique 9 
together with Eqs. (39) and (40), as listed in Table I. 
For comparison, Table I also lists the first line of each 
spectroscopic series, i.e., the circular lines from quan- 
tum level level k + 1 to quantum level k. Historically, 
the series of hydrogen were named after Lyman, Balmer, 
Ritz-Paschen, Brackett, etc. The frequency over reduced 
mass of the first line of each spectroscopic series in atomic 
units is the second column of Table I. We used a Newton 
method in the complex-A plane to solve Eqs. (39), (40) 
and (55), as used in Ref. [15] with Dirac's theory. Our 
calculations for hydrogen used the protonic-to-electronic 
mass ratio (m 2 /m 1 ) = 1836.1526. Table I gives fre- 
quency over reduced mass calculated by QM for the first 
line of the k th spectroscopic series (atomic units), numer- 
ically calculated orbital frequency over reduced mass (for 
a suitable integer q(k)), (137 3 0)/^ = 137 3 9 2 (e xy - e z ), 
angular momentum of unperturbed circular guide, and 
integer q. 

Table I is to be compared with Table I of Ref. 
[15], which discusses Dirac's electrodynamics with self- 
interaction. The same surprising agreement is found in 
Ref. [15], even though Dirac's theory has limitations not 
shared with the Wheeler-Feynman theory, as discussed 
in Ref. [2]. Notice we also jumped q = 10 in Ref. [15], 
again because the resonance condition is only necessary. 
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TABLE I. Numerical calculations for hydrogen with 
(m2/mi) = 1836.1526. Quantum number k of the circular 
Bohr transition k + 1 — > k, frequency over reduced mass of 
the circular QM line in atomic units, wqm = — (k+i)' 2 )' 
orbital frequency over reduced mass in atomic units, = 
137 3 6> 2 (e xy — angular momentum of unperturbed orbit in 
units of e 2 /c, l z = and integer q. 

As mentioned below Eqs. (48) and (49), the ener- 
getic mismatches are monotonically increasing for q > 6, 
and Table I frequencies are within ten percent agree- 
ment with the first 12 circular hydrogen lines starting 
from q = 7 at k = 1, i.e., q(l) — 7. Since condition 
(55) is only necessary (and not sufficient), some values 
of q may correspond to unstable orbits. This is analo- 
gous to the description by QM [9], where there are se- 
lection rules on top of three conditions involving integer 
quantum numbers, and Table II includes the lines that 
were skipped in Table I. A theory for the g's that were 
skipped is presently lacking. Inspection of Table II shows 
that for q — 1,2, 3, 4, 5, 6 condition (55) predicts 9 still in 
the atomic range, but the angular momentum spacing is 
about half of Planck's constant. Table II also includes the 
skipped numerical calculations for q = 8 and q = 10. A 
suggestion from Sommerfeld's quantization is that there 
must be three conditions like (55), involving three inte- 
ger quantum indices, not just one, reinforcing that our 
single condition (55) is only necessary and alone might 
not determine a stable orbit. Inspection shows that if 
the values of Table II were included in Table I, the an- 
gular momentum jump from consecutive lines would be 
much lesser than about a hundred units of e 2 /c, which is 
suggestive of what the missing conditions should do. 

Next we show the numerical calculations for the muo- 
nium atom using the positive-muon-to-electron mass ra- 
tio (mj/mi) = 1836.1526/9. Table III lists the frequency 
over reduced mass of the first line of each spectroscopic 
series as calculated with QM (in atomic units), orbital 
frequency in atomic units, (137 3 f2)/^ = 137 3 9 2 (e xy ~ e z ), 
and angular momentum of unperturbed circular orbit in 
units of e 2 /c. The agreement of the numerical calcula- 
tions with nature and QM is again within a few percent 
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TABLE II. Numerical calculations for hydrogen with q < 7 
and (m2/mi) = 1836.1526. orbital frequency over reduced 
mass in atomic units, (fi//- 1 ) = 137 3 # 2 (e xy — e z ), angular 
momentum of unperturbed circular guide in units of e 2 /c, 
lz = 0' 1 , and integer q. 



for frequencies. 



k 


WQM 


lZ7 6 6 z {e xy -e z ) 


lz = e~ L 


q 


1 


3.750xl0 _i 


4.039X10 -1 


185.63 


7 


2 


6.944xl0~ 2 


8.762xl0~ 2 


308.94 


9 


3 


2.430 xlO" 2 


2.356xl0~ 2 


478.65 


11 


4 


1.125xl0" 2 


1.304xl0" 2 


582.96 


12 


5 


6.111xl0~ 3 


7.491xl0~ 3 


701.31 


13 


6 


3.685xl0 -3 


4.445 xl0~ 3 


834.53 


14 


7 


2.406 xlO" 3 


2.717xl0 -3 


983.41 


15 


8 


1.640xl0~ 3 


1.704xl0~ 3 


1148.75 


16 


9 


1.173xl0 -3 


1.095xl0~ 3 


1331.34 


17 


10 


8.678 xlO -4 


7.187xl0 -4 


1531.19 


18 


11 


6.600xl0~ 4 


4.810xl0~ 4 


1751.36 


19 


12 


5.136xl0~ 4 


3.277xl0~ 4 


1990.34 


20 



TABLE III. Numerical calculations for muonium with 
(7712/mi) = 1836.1526/9. Quantum number k of the circu- 
lar Bohr transition k + 1 — > k, frequency over reduced mass of 
the circular QM line in atomic units, wqm = \ — (k+i)' 2 )' 
orbital frequency over reduced mass in atomic units, (£1 / p.) = 
137 s 8 2 (e xy — Ez), angular momentum of unperturbed circular 
guide in units of e 2 /c, l z = 0~ 1 , and integer q. 

We repeat the numerical calculations for positronium 
using the positron-to-electron mass ratio (mz/mi) = 1: 
Table IV lists the frequency over reduced mass of the 
first line of each spectroscopic series calculated by QM 
(in atomic units) , orbital frequency over reduced mass in 
atomic units and the angular momentum of the unper- 
turbed circular orbit in units of e 2 /c. Notice in Table IV 
that for positronium the values of l z = 1/9 are consis- 
tently larger. Using again the fact that condition (55) is 
only necessary, Table IV starts the q when the angular 
momentum spacing is about constant, i.e., q — 8. For 
positronium the numerical calculations find the first root 
1/9 = 40.501 only at q = 3, again in the atomic magni- 
tude. The spectrum agrees with QM within less than a 
few percent for the circular lines of the first 12 series. 

In the days of Bohr, only about 12 lines of the Balmer 
series could be observed with vacuum tubes, and about 
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TABLE IV. Numerical calculations for positronium with 
(7712/mi) = 1. Quantum number k of the circular Bohr 
transition k + 1 — ► fc, frequency over reduced mass of the 
circular QM line in atomic units, wqm = §(^i — (gjlp )i or ~ 
bital frequency over reduced mass in atomic units, (f2//x) = 
137 3 6> 2 (e aiy — £z), angular momentum of unperturbed orbit in 
units of e 2 /c, l z = 6 1-1 , and integer 5. 



for q = 7 and (jj,/M) = 9/1836 (muonium), Eq. (56) 
yields 9 ~ 187.17, and last for g = 7 and (ji/M) = 1/4 
(positronium), Eq. (56) yields ~ 186.81. The numer- 
ically calculated values of l z = 6~ x fall approximately 
between the two consecutive Bohr orbits k and k + 1 (ar- 
tificially) used to calculate frequencies, and Tables I, III 
and IV show that frequencies of spectral lines agree even 
better with QM and nature. 

Last, the orbital frequencies determined by the 
Weierstrass-Erdmann mechanism (55) can be expressed 
as a difference of two spectroscopic terms, just like in 
QM, i.e., 



3 6e xy {q) -6e z (s) 



(57) 



with spectroscopic terms Oe xy (q) and 6e z (s) defined by 
the eigenvalues of two linear and infinite dimensional 
problems, i.e., Eq. (31) and Eq. (B16) of Appendix B in 
Ref. [15]. 



thirty-three from spectra of celestial bodies [26]. Sur- 
prisingly, the emission frequencies agree better with QM 
for the decays from the 12 deepest quantum levels. As 
explained in [15], the cancelation of dipolar far-fields 
involves quadratic terms that might require larger am- 
plitudes at large q. Given that the xy modes modify 
the unperturbed z-angular momentum more and more 
at larger q, our perturbative results should get worse 
at larger q. We carried our numerical calculations up 
to q = 43 (not shown in Table II), and observed that 
the angular momentum spacing increases slowly with 
q. By defining an effective k = (Z z /170), the agree- 
ment of the numerical calculations with an effective angu- 
lar momentum separation seems to continue indefinitely 
within thirty percent, suggesting we can approximate a 
large number of eigenvalues near the discrete spectrum 
of Schroedinger's equation. The agreement with emission 
lines above k = 13 slowly deteriorates, which is suggestive 
that the corresponding broken minimizers are becoming 
far from planar. A one-to-one comparison with natural 
spectra should wait the investigation of broken extrema 
with spikes filling a tridimensional region. 

The surprising agreement of the numerical calculations 
with a universal value for the fine-structure co nstan t is 
due to the logarithmic dependence of a = ln(y^pO) a s 
explained above Eq. (37). Moreover, because of the weak 
dependence on 9, Eq. (56) is an implicit transcendental 
condition and the 0(1/ X 2 ) formulas (48) and (49) yield 



1 



1 



(a 2 + 7rV ) 2 
8iraq 



(56) 



The roots of (56) are insensitive to changes in (/i/M) over 
three orders of magnitude, e.g., for q — 7 and (/x/M) = 
1/1837 (hydrogen), Eq. (56) yields 6 ~ 190.09, while 



GAC AND VANISHING FAR-FIELDS 

Radiation in Maxwell's electrodynamics is treated with 
Poynting's theorem [5], which within variational electro- 
dynamics holds only for extended fields of globally C 2 
orbits, in which case Maxwell's electrodynamics is a spe- 
cial case of variational electrodynamics. In the follow- 
ing we use only generally valid properties of variational 
electrodynamics. The physical meaning of the extended 
far-fields in variational electrodynamics is only that of a 
coupling to distant charges. The significance of a surface 
integral of the Poynting vector of extended fields unde- 
fined in sets of zero measure is an open application of the 
divergence theorem of modern functional analysis [10]. 

Natural spectra of atomic gases containing Avogadro's 
number of atoms pose a many body problem where each 
atom suffers perturbations from all other atoms. This un- 
less a decoupling from other atoms or laboratory bound- 
aries is achieved by the GAC. In Ref. [2] it was shown 
that globally bounded two-body orbits with vanishing 
extended far-fields must involve velocity discontinuities 
along the continuous trajectories. 

Under the above physical interpretation of extended 
far-fields as couplings, we recall the far-field- vanishing 
mechanism of Ref. [15] using the quadratic correction to 
far-field (4) due to the last right-hand-side term of (5), 



J J f ., 



(58) 



Not surprisingly, quadratic term (58) is ignored when 
approximating the two-body motion with a Coulombian 
ODE. Again, Colombian ODE's have a tangent dynam- 
ics spanned by a finite number of slow frequencies [25], 
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which cause small accelerations, while the tangent dy- 
namics defined by NDDE (31) has arbitrarily large fre- 
quencies and arbitrarily large accelerations. Since ac- 
celerations are arbitrarily large and transversal veloci- 
ties quasi-luminal, the quadratic term (58) introduces 
beats of the fast frequencies [15]. The GAC mechanism 
of Ref. [15] to cancel the dipolar far-fields of the unper- 
turbed (slow) circular motion by (58) with a resonant 
beat among fast planar and fast perpendicular modes 
would be inexistent for a Coulombian ODE. 

It is a nice surprise that the WECC are compatible 
with the GAC. In fact, (55) with q = s is equivalent to 
the condition for the quadratic term to cancel the dipolar 
far-fields of the unperturbed circular orbit [15]. Equa- 
tion (55) with q — s is exactly Eq. (56) of Ref. [15]. 
Therefore, (55) determines the finitely measured family 
of orbits that can acquire vanishing far-fields by linear 
perturbations, as necessary to decouple from experimen- 
tal boundaries or other atoms of the gas. 

The logical application of multiscale perturbation is 
to balance the fast timescale first. For that one could 
freeze the slow guiding motion more generally than as 
a circular orbit. As long as the slow velocity is of the 
order determined by the centripetal force of the Coulom- 
bian limit, i.e., [iv 1 ~ 1/r, the far-field mechanism of 
0(1/A 2 ) corrections to (33) and (36) prescribes the same 
resonance condition (55). This is further indication that 
our determination of atomic scales is not accidental, and 
other finitely measured families of orbits with vanishing 
far-fields should be found again in the atomic magnitude. 

A second confirmation that the WECC mechanism se- 
lects the atomic scale is that the simple 180-degree re- 
bouncing mechanism turns off at smaller scales. Ac- 
cording to Eq. (51), 180-degree corners are an alterna- 
tive to opposite charges falling into each other only if 
f > , 1 at which distance the mechanism ceases to 

— v /mim 2 ' 

operate, thus creating a confinement core at the nuclear 
maqnitude f < , 1 

We do not claim to have found all possible spiky orbits, 
and the q = s family of periodic orbits with vanishing far- 
fields precisely at the atomic magnitude suggests a sys- 
tematic search should continue. For non-harmonic (e.g. 
elliptical) orbits, an effective mechanism of cancelation 
of the main harmonic of the slow periodic guide should 
operate for q ^ s. Investigations of stability are better 
done numerically, and for that finite element methods are 
being developed [18, 23]. 

DISCUSSIONS AND CONCLUSION 

Here, as in [18], we adopted particle time for each 
particle's evolution parameter. The Fokker action 
with two arbitrary evolution parameters [1], is im- 
proper for the finite-element method of [18] because 
the parametrization-invariance symmetry causes a zero 



eigenvalue and a singular Hessian matrix [18]. WECC 
from the parametrization-invariant action were used 
in [27] to model double-slit interference caused by 
interaction-at-a-distance with the corner points of the 
trajectories of material electrons inside the grating [27]. 
The parametrization-invariant action yields four scalar 
conditions for the WECC1, while WECC2 is trivially sat- 
isfied [28] . The choice of particle-time evolution param- 
eter yields three scalar conditions for the WECC1 and 
one scalar condition for WECC2 (an equivalent set of 
four conditions as explained in [28]). 

The GAC uncoupling mechanism operates even at 
closer distances along globally bounded trajectories of 
atoms with a zero net charge, because then also the near 
fields vanish as a consequence of Gauss's integral law. 
Along trajectories satisfying the GAC, atoms with a zero 
net charge can stand decoupled at distances ranging from 
intermediate to far distances. The minimizer for the com- 
bined system is the direct product of the unperturbed 
trajectories precisely because each atom produces no far 
fields and little near fields to act back on the other atoms. 

A motivation for Wheeler-Feynman theory [6] was 
Sommerfeld's quantization conditions [20, 21, 29]. Sur- 
prisingly, broken extrema of variational electrodynamics 
satisfy a condition involving precisely the same relevant 
quantities of quantum mechanics, i.e., momenta and en- 
ergies [4]. Energy and momentum entered Somerfeld's 
theory via Hamiltonian mechanics, which became para- 
doxal after the no-interaction theorem proved that only 
trivial non-interactive two-body motion is covariantly de- 
scribed by a Hamiltonian [16, 30, 31]. 

Variational electrodynamics is also free of the 
Aharonov-Bohm paradox [32] along broken extrema. The 
origin of this paradox is that along C°° orbits the elec- 
tromagnetic equations of motion involve only derivatives 
of the vector potential: The Euler-Lagrange equation of 
partial Lagrangian (15) yields (1) with a right-hand-side 
equal to VC/j + dAj/dt - x (V x A-j)- Instead, along 
broken extrema the vector potential Aj itself appears on 
the WECC1, Eq. (14), determining an interference-at-a- 
distance just like in QM, as discussed in [27]. 

There is no fitting parameter in our calculations, and 
the agreement with nature proves variational electrody- 
namics is a physically sound extension. Agreement be- 
yond thirty percent should wait a perturbative scheme 
with corrections for the guiding center motion. 

The general hydrogenoid spectrum falls in the class 
Hqm = 2 (p- — js) [9]. If variational electrodynamics is 
to approximate the whole hydrogenoid spectrum, broken 
minimizers still to be discovered should correspond to the 
elliptical lines (j ^ k + 1). For example, spiky orbits with 
two guiding circles, with each particle diving to a lower 
circle and riding along it for a while, then jumping back 
to an upper circle and riding along it for a while, and so 
on. These orbits involve two Weierstrass-Erdmann con- 
ditions and two integers, just like the hydrogenoid series 
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[9, 26]. A classification of spiky orbits and a multiscale 
scheme starting from fast motions about non-smooth and 
almost everywhere slow guiding motions is lacking. 
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